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Hyperbolic  Systems  of  Conservation  Laws  in  Several  Space  Variables 


Peter  D.  Lax 

Courant  Institute  of  Mathematical  Sciences,  New  York  University 


Introduction 

In  Section  1  we  show  how  to  associate  to  each  analytic 
function  f  of  one  variable  a  complex  conservation  law.  When 
written  in  terras  of  real  variables,  this  equation  turns  into  a 
symmetric  hyperbolic  system  of  two  first  order  conservation  laws 
in  time  and  two  space  variables  for  two  unknowns.  The  hope  is 
that  the  theory  of  analytic  functions  can  be  exploited  to  study 
the  properties  of  the  solutions  of  these  model  equations. 

These  systems  are  strictcly  hyperbolic  where  f  ^  0,  and 
nonlinear  where  f"  ^  0.  There  is,  however,  in  every  state  at 
least  one  direction  of  propagation  which  is  not  genuinely 
nonlinear. 

In  Section  2  we  show  that  every  2x2  quasilinear  first 
order  hyperbolic  system  fails  to  be  genuinely  nonlinear  in  some 
direction  of  propagation. 
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The  notions  of  the  first  section  are  asily  extended  to 
analytic  functions  of  several  variables. 

This  paper  is  dedicated  to  Sigeru  Mizohata,  in  recognition 
of  his  deep  contributions  to  our  understanding  of  the  theory  of 
partial  differential  equations. 

1 .    Symmetric  Hyperbolic  Systems  of  Conservation  Laws. 

f  =  f(U)  denotes  an  analytic  function  of  a  single  complex 
variable  U.  We  impose  on  the  complex  valued  function  U(z,t)  of 
the  complex  variable  z  and  the  real  variable  t  the  following 
nonlinear  partial  differential  equation: 

(1°1)  S^U  +  3zf(U)  =  0  . 


Here,  as  customary,  9^  denotes 


(1-2)  ^z   =  J  (\   -  ^    9y).     z  =  X  +  iy 

and   the   bar  denotes  the  complex  conjugate.   We  can  express  this 
equation  in  terms  of  the  real  and  imaginary  parts  of  U  and  f,,: 


(1-3)  U  =  u  +  iv  ,  i  fy  =  a  +  ib 


Setting  these  into  (1.1)  gives 
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(1.4) 


u,.  -  iV(.  +  (a+ib)(aj^+iv^-iu  +v  )  =  0  , 


The  matrix  form  of  this  equation  is 


(1-^)'  Ot   ^^O^^^lVy   =  0 


where 


(1.5) 


A  = 


a  -b' 
-b   -a 


b   a 

a   -b 


Since  A  and  B  are   symmetric   matrices,   (1.4)'   is   a   symmetric 
hyperbolic  system. 

The   characteristics   of   this   sytem  are  the  eigenvalues  of 
5A  +  nB;  these  are 


(1.6) 


/a^+b^   for   ?2+^2  ^  ^ 


Thus  the  speed  with  which   signals   propagate   from  a   point   is 


^a   +h      in  all  directions. 

It  was  observed  in  [3],  see  also  [5]  and  [2],  that  a 
symmetric  system  of  conservation  laws  satisfies  an  additional 
conservation  law.  This  can  be  derived  from  the  complex  form 
(1.1)  of  the  equation  by  multiplying  it  by  U: 
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(1.7)  UU^.  .+  UfyU^  =  0  . 

Denote  by  h   the   analytic   function,   determined   up   to   a 
constant,  whose  derivative  is  f: 


(1.8)  hy   =  f 


Then 


(1.8)'  Ufy  =   (Uf   -  h)u  ■-     ky 


Equation  (1.7)  can  be  rewritten  as 


(1.8)"  UU^  +1^2  =  0 


Add  the  complex  conjugate  of  (1.8)  to  (1.8"): 


(1.9)  MUI^  +  k^  +  k_  =  0 

y 


|U|  is  a  convex  function  of  U;  as  explained  in  [5],  it  plays  the 
role  of  an  entropy  for  equation  (l.I),  in  the  following  sense:  if 
U  is  a  strong  limit  as  e  -»■  0  of  solutions  U^^^  of  the  viscous 
equations 
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(1.10)  U(^)  +  f(u(0)^  =  tAU^^l 


then  U  satisfies  the  entropy  inequality 

(1.11)  3^|U|2+lc^+k~<0. 

nil      z     -  _ 

At  a  point  of  discontinuity  of  U,  called  a  shock,  the  strict 
inequality  holds  in  (1.11). 

We  note  that  every  antianalytic  function  of  U  is  a  tirae- 
independent  solutiion  of  (1.1). 

We  show  now  how  to  associate  to  each  analytic  function  h  of 
several  complex  variables  a  symmetric  system  of  conservation 
laws . 

Let  h(U  ,...,U'^)  be  an  analytic  function;  denote  its  partial 
complex  derivatives  by 


(1.12)  fJ  =  h  . 


The  system  of  n  complex  conservation  laws 


(1.13)  UJ  +  fj^  =  0  ,   j  =  l,...,n, 

when  written  in  terms  of  the  real  and  imaginary  parts  of  L'J  ,  is  a 
symmetric  system  of  conservation  laws;  this   is   easy   to   verify 
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directly,  but  it  is  still  easier  to  verify  an  equivalent  fact, 
that  Z  I'j'-'l  satisfies  an  additional  conservation  law.  To  see 
this  multiply  the  j^"  equation  (1.13)  by  UJ  and  sum  over  j: 

(1.14)  E  UJ  UJ  +  E  UJf J  =  0 

Using  (1.12)  we  see  that  the  second  term  can  be  written  as 


(1.15)  k^  =  (Z  UJfJ  -  h)^  . 


Taking  the  real  part  of  (1.14)  and  using  (1.15)  gives 


(1.16)  8^  E  |UJ  |2  +  k^  +  k_  =  0 


z 


Equation   (1.16)   holds   for   all  smooth  solutions  of  (1.13);  for 
discontinuous  solutions  the  corresponding  inequality 

(1.16)  3^  E  luJ |2  +  k^  +  k7  =  0 

z 

serves  as  entropy  condition. 

2.    Linear  Degenerations  of  Quasilinear  Equations 

A  first  order  system  in  t  and  one  space  variable, 

(2.1)  u^  +  Au^  =  0  , 
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a  a  vector  function  with  values  in  1",  A  an  nxn  real  matrix,  is 
called  strictly  hyperbolic  if  the  matrix  A  has  n  real  and 
distinct  eigenvalues  X^,...,x^.  When  (2.1)  is  quastlinear,  i.e. 
when  A  depends  on  u,  so  do  the  eigenvalues  \,,...,X  and  the 
corresponding  right  eigenvectors  r,,..,,r  : 

(2.2)  Ar  =  Xr  . 

The  system  (2.1)  is  called  genuinely  nonlinear  if   for  each 
eigenvalue   A^   =   x,   the   gradient  of  X  satisfies  for  all  u  the 

inequality 


(2.3)  A^  .  r  ?t  0  , 


where  r  is  the  corresponding  right  eigenvector.   Note  that   (2.3) 
is  independent  of  the  normalization  r. 

For  n  =  1,  (2.1)  is  a  scalar  equation: 

u^  +  a  Uj^  =  0  , 

a   =  a(u)   =  A(u).   In  this  case  there  is  a  single  eigenvalue,  a 
itself,  and  condition  (2.3)  requires  that 


(2-3)'  -,^0, 
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clearly  the  condition  of  genuine  nonlinearity . 

The  significance  of  genuine  nonlinearity  for  the  theory  of 
hyperbolic  equations  and  the  formation  of  singularities  is- 
explained  in  any  treatise  on  the  subject  such  as  [6]  or  [8]. 

We  turn  now  to  first  order  systems  in  two  space  variables: 

(2.4)  u^.  +  Au^  +  Bu   =  0  , 

u  an  1"  valued  function  of  x,y,t,  A  and  B  real  n  x  n  matrices. 

Equation  (2.4)  has  solutions  u  that  propagate  in  the 
direction  (5>n),  i.e.   are  of  the  form 

(2.5)  u(x,y,t)  =  w(Cx  +  ny,  t) 
Such  a  u  satisfies  (2.4)  if  w(s,t)  satisfies 

(2.6)  w^  +  C  Wg  =  0  , 
where 

(2.7)  C  =  ?A  +  nB  . 

Equations  (2.4)  is  called  strictly  hyperbolic  if  equation 
(2.6)  characterizing  those  solutions  that  propagate  in  the 
direction   (C,n)   is   strictly  hyperbolic  for  every  choice  of  the 
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direction.  That  means  that  the  matrix  C  has  to  have  real  and 
distinct  eigenvalues  for  all  real  choices  of  E,,    n,  5  +n  >  0. 

Suppose  equation  (2.4)  is  quasilinear,  i.e.,  that  A  and  B 
are  functions  of  u.  We  call  (2.4)  genuinely  nonlinear  if  the 
equations  (2.6)  for  unidirectional  waves  are  genuinely  nonlinear 
for  every  choice  of  the  direction. 

Take  a  scalar  equation 

^t     *  ^"X     +     t)Uy     =     0      , 

a,b  functions  of  u;  equation  (2.6)  is 

W^_  +  C  Wg  =  0  , 

c  =  Ca  +  nb  . 
Condition  (2.3)'  becomes 


Ca^  +  nb^  ?t  0  , 


clearly  violated  by  suitable  choice  of  (C,'^).  Thus  a  scalar 
quasilinear  equation  in  two  space  variables  is  never  genuinely 
nonlinear  in  all  directions. 

We  turn  now  to  the  2x2  system  (1.4)'.   We  write 


i.e. 
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a  +  ib  =  qe^'t' 

a  =  q  cos  <i>    ,      b  =  q  sin  <^    , 


so  that  the  matrices  A,  B  in  (1.5)  can  be  written  as 


A  =  q 
Setting 


cos  <p ,    -sin 
-sin  (Ji ,  -cos 


B  =  q 


sin  4),   cos  4) 
cos  ■}) ,   -sin  (p 


E,   =  cos  6  ,   n  =  sin 


we  get  for  C  given  by  (2.7) 


D  =  q 


cos  (e-(|))   sin(9-^) 
sin  (B-'p)    -cos(9-(t)) 


The  eigenvalues  of  C  are  ±  q;  a  brief  calculation  shows  that   the 
right  eigenvectors  are 


(2.8)       r^  = 


cos 


sm 


9-(j) 


9-^ 


,   r. 


sm 


-cos 


9-^ 


Condition  (2.3)  in  this  case  is 
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We  claim  that  this  cannot  hold  for  all  6.  For  q,  and  therefore 
q^j  is  independent  of  9;  on  the  other  hand  formulas  (2.8)  show 
that  as  9  goes  from  0  to  2ti,  the  vectors  r^  and  r_  turn  over  an 
angle  of  tt  ;  therefore  there  are  values  9^.  when  r+(9+)  are 
orthogonal  to  g^  .  For  these  values  of  9  condition  (2.9)  is 
violated. 

Andy  Majda  raised  the  question  if  there  genuinely  are 
nonlinear  2x2  systems  in  two  space  dimensions.  It  turns  out 
that  indeed  there  aren't  any: 

Theorem.  Every  real,  strictly  hyperbolic  quasilinear  system 
(2.4)  for  two  unknowns  of  two  space  variables  is  linearly 
degenerate  in  some  direction. 

Proof :  Define 

(2.10)  C(9)  =  cos9  A  +  sine  B  . 

Denote  the  eigenvalues  of  C(8)  by  X+{Q),    labeled  so  that 

(2.11)  X_(9)  <  X^(9)  ; 
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the  corresponding  right  eigenvectors  are  denoted   by   r_^(9).    We 
choose  them  to  be  real  and  of  unit  length: 

(2.12)  |r+(e)|   =  1  , 

in  any   old   norm.    This   still  leaves  an  arbitrary  factor  ±  1, 
which  we  fix  arbitrarily  at  6  =  0,   and   for  all  other   6  in 
0  <  6  <  2tt  by  requiring  r  +  (9)  to  vary  continuously  with  9. 
The  matrix  C(9)  defined  by  (2.10)  satisfies 

(2.13)  c(e  +  IT)  =  -  c(e)  . 

It  follows  from  (2.13)  and  (2.11)  that 

^+(9  +  tt)  =  -  X_(9) 
X_(9  +  IT)  =  -  \  +  (9)  . 

It  follows  from  this  and  (2.12)  that 

r^(9+Tr)  =  a^r_(9) 
(2.14) 

r_(e+it)  =  a_r+(9) 

where   a_j_  and   a_  are   ±  1.    Since   r+(9)   were   chosen   to  be 

continuous  functions  of  9,   a+  also   are   continuous   functions; 

since  their  value  is  ±  1 ,  they  are  constant. 
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Since   r^.(9)   vary   continuously,   the   orientation  of   the 
ordered  base 

(r_(8),  r^(9)} 
does  not  change;  in  particular 

{r_(0),r  +  (0)}   and   {r_(  ir) ,  r+(  ti)  } 


have  the  same  orientation.  By  (2.14) 


fr_(0),r+(0)}   and   {a_r+(0) , a^r_(0)} 


have  the  same  orientation;  this  proves  that 


(2.15)  a^o.  =  -  1 


Setting  9  =  0  and  9  =  tt  we  obtain 


r+(2TT)  =  a^.r_(TT)  =  a+a_r  +  (0)  ; 


so  by  (2.  15) 


(2.16)^  r^(2TT)  =  -  r^(0) 
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and  similarly 

(2.16)  r_(2TT)  =  -  r_(0)  . 

The  eigenvalues  X+(u,9)  are  periodic  functions  of  6  with   period 
2iT ;  therefore  so  is  their  gradient;  this  and  (2.16)  shows  that 


(2.17)  ^u(2^)  •  '^(^rr)  =  -  X^(0)-r(0) 


for  both  eigenvalue-vector  pairs.  Since 

(2.18)  X^(9)  .  r(e) 

varies  continuously  with  9,  it  follows  from  (2.17)  that  (2.18) 
vanishes  for  some  9.  By  (2.3)  the  system  (2.4)  fails  to  be 
genuinely  nonlinear  in  the  direction  (cos  9,  sin  9).  o 

The  argument  presented  above  shows  this 

Corollary.  Let  (2.4)  be  a  2n  x  2n  real,  quasilinear, 
strictly  hyperbolic  system  in  two  space  variables,  n  odd.  Then 
(2.4)  fails  to  be  genuinely  nonlinear  in  some  direction. 

The  topological  argument  used  in  the  above  proof  is  taken 
from  [7],  where  it  is  used  to  prove  that  there  are  rm  strictly 
hyperbolic  systems  of  2n  equations  for  2n  unknowns  in  three  space 
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vartables  for  n  odd.   The  interesting  question  of   the   existence 
of   m  X  m  strictly   hyperbolic   systems   in   three   variables  is 
handled  for  arbitrary  m  in  the  interesting  paper  [1]. 
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